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We consider the dynamics of electromagnetic fields in an almost - Friedmann- Robertson - Walker 
universe using the covariant and gauge - invariant approach of Ellis and Bruni. Focusing on the 
situation where deviations from the background model are generated by tensor perturbations only, 
we demonstrate that the coupling between gravitational waves and a weak magnetic test field can 
generate electromagnetic waves. We show that this coupling leads to an initial pulse of electro- 
magnetic waves whose width and amplitude is determined by the wavelengths of the magnetic field 
and gravitational waves. A number of implications for cosmology are discussed, in particular we 
calculate an upper bound of the magnitude of this effect using limits on the quadrapole anisotropy 
of the Cosmic Microwave Background. 
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I. INTRODUCTION 

There have been numerous investigations on the scat- 
tering of electromagnetic waves off gravitational fields 
(See Refs. jl| for a representative sample). Most of this 
research have been focused on the effect gravitational 
waves have on vacuum electromagnetic fields. Other pa- 
pers, see e.g. Refs. ||, also consider situations in Astro- 
physics where plasma effects are taken into account. 

Magnetic fields play an important role in our Universe, 
appearing on all scales from the solar system, through 
interstellar and extra - galactic scales, to intra - cluster 
scales of several Mpc. Although magnetic field inhomo- 
geneities have not yet been observed on scales as large 
as those exhibited by Cosmic Microwave Background 
(CMB) anisotropies, it is natural to expect that mag- 
netic fields exits on such scales ||, and that they could 
play a role in the formation of large-scale structure. In- 
deed many mechanisms have been proposed to explain 
how these fields may be generated in the early universe - 
a process called primordial magnetogenesis. For example 
on small scales (less than the Hubble radius), QCD and 
Electro -Weak transitions can give rise to local charge 
separation leading to local currents which can generate 
magnetic fields |Q . Large - scale magnetic fields can be 
generated during inflation or in pre -Big Bang models 
based on string theory |§, in which vacuum fluctuations 
are amplified via the inflaton or dilaton. 

The effect magnetic fields have on density perturba- 
tions has been studied extensively by a number of au- 
thors, both in the context of Newtonian and Relativistic 
Cosmology |||7j , but as yet there have been no studies of 
their effect on gravitational wave perturbations. 

In what follows, we use the well known covariant and 
gauge - invariant approach of Ellis and Bruni M to in- 



vestigate this interaction in the context of cosmology by 
considering the dynamics of electromagnetic fields in an 
almost Friedmann - Robertson - Walker (FRW) universe, 
focusing in the situation where deviations from the FRW 
background are generated by tensor or gravitational wave 
perturbations P,|i~o||. 

We show that in the presence of a weak (near- 
homogeneous) magnetic test field propagating on the 
background FRW model pl[ 3 the gravitational waves 
couple non- linearly to this field to produce a pulse of 
gravitationally induced electromagnetic waves. In par- 
ticular, because of the different ways in which tensor per- 
turbations enter the wave equations for the electric and 
magnetic fields, respectively, there will be, in the case 
of long wavelength gravitational waves and large - scale 
magnetic fields, a growth in the expansion normalised 
electric field, as the expansion normalised shear grows in 
time. 

This paper is organised as follows. After a short dis- 
cussion of notation and conventions, in section [II wc 



outline in detail the linearisation procedure used to ap- 
proximate the Einstein - Maxwell equations in cosmol- 
ogy. In section [V we derive a set of non -linear wave 



equations which show how electromagnetic fields can be 
generated when gravitational waves couple to an near- 
homogeneous magnetic test field propagating on the 
FRW background. Finally in section |y| we solve the equa- 
tions perturbatively and use our results to put an upper 
bound on the size of these gravitationally induced fields. 



II. NOTATION AND CONVENTIONS 

Notation and conventions are taken to be the same 
as in p2]. In particular 8itG = c = 1; the projected 
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spatial covariant derivative a tensor T cd " e /... is given 
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where u a is the 4 - velocity of the matter, g ab is the metric 
tensor and h ab = g ab + u a u b is the spatial projection ten- 
sor (h ab u b = 0). A dot denotes the covariant derivative 
along u a , so for any tensor T cd " e /... = u a V a T cd - " e /...- 
We assume that the matter is described by irrotational 
dust | [r^ | so that the pressure p, acceleration vector u a and 
vorticity tensor uj a b all vanish exactly. In this case the 
first covariant derivative of the 4 - velocity can be written 
as V Ufc = cr ab + ^Qh ab , where a ab and 9 are respectively 
the usual shear and volume expansion of the matter con- 
gruence. We also define the Hubble parameter H terms 
of the expansion 8 and scale factor a in the usual way: 
H = 6/3 = a/a. 



like the tensor perturbations of the Bardeen approach 
0- 

Given the assumed equation of state, these conditions 
also imply that the spatial gradient of the expansion V a 
vanish (see ||). Together with these conditions pro- 
vide a unique characterisation of tensor perturbations. 



IV. EINSTEIN - MAXWELL EQUATIONS 

We assume overall charge neutrality and use the 
Bianchi identities and Maxwell's equations as presented 
in Jl^]. Then, with the above prerequisites, we obtain 
a set of non - linear wave equations for the gravitational 
(a ab : shear) and electromagnetic (E a , B a : electric and 
magnetic fields) degrees of freedom: 



III. APPROXIMATIONS 

In order to simplify the non -linear dynamics of the 
coupled Einstein - Maxwell equations and to isolate the 
effects we are looking for, we will adopt the following 
approximation scheme based on two parameters: e g will 
refer to quantities occurring in the gravitational equa- 
tions, while £ e m characterise the electromagnetic field. 
We assume that the gravitational equations follow the 
almost FRW conditions |8), so that the energy density 
/i and expansion have a non - zero contribution in the 
background model and can therefore be considered O(0 g ) 
while a a b, E ab , and H ab vanish in the background and 
are 0(e g ). In the case of the Maxwell field, we assume 
that there is a weak magnetic test field Bq at O(0 om ) 
which propagates on the background FRW model, whose 
gravitational influence is given by the Alfven parameter 
e = (B^Bl/ii) 1 / 2 . On the other hand the electric field 
E a vanishes in the background and is considered to be 
0(e em ). The perturbation scheme we adopt is to drop 
terms of 0(e 2 ) (so that the magnetic field does not con- 
tribute to the gravitational dynamics [Q), 0(e 2 ), 0(s 2 m ), 
and C(e g £cm)- 

In the covariant approach to linear perturbations of 
FRW models pure tensor or gravitational wave per- 
turbations are characterised by the following covariant 
conditions |g|JTc||: 



V b E ab = 
V b H ab = 



=> V o m = 0, 

=> U a =0., 
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the first one excludes scalar (density) perturbations and 
the second, vector (rotational) perturbations. The con- 
ditions that the terms on the right hand side vanish, are 
analogous to the transverse condition on tensor pertur- 
bations in the metric approach. In addition, we notice 
that since the Weyl tensor is the trace - free part of the 
Riemann tensor, both Eat, and H ab are trace -free, again 
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where 
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are gravitational induced magnetic and electric currents, 
and E ab = -& ab - f Qa ab and H ab = e cd( - a V c a% are 
the electric and magnetic parts of the Weyl tensor. Also 
A/ = / — V 2 / where / is any tensor orthogonal to u a . 

In the above equations, the electric and magnetic fields 
consist of two parts, a contribution due to the magnetic 
test field Bq which gives rise to the current j a [|l5| and 
contributions generated by the non -linear coupling of 
this test field to gravitational waves via the gravitation- 
ally induced currents j% and j% : 



E a = EL. 



B a = Bq 
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V. ANALYTIC SOLUTIONS AND NUMERICAL 
INTEGRATION 

We solve the above equations perturbatively by first 
calculating the gravitationally induced currents j% and 
j% and then solving (2b) and (^) together with (pa|) for 
the gravitationally induced electric and magnetic fields. 

To O(0om), Maxwell's equations @ give E — Q and 
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which we can integrate to obtain: 
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where A? » and .A"^ determine the spatial variation of 
the magnetic test field and are constant along the fluid 
flow lines: A a {n) = A a A =00. Furthermore we as- 
sume that the spatial functions AJ, and A^K satisfy the 
Helmholtz equation 
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in this way defining a specific length scale Xb — 2-n:a/n 
associated with the magnetic field determining its scale of 
inhomogeneity, where n is a fixed wavenumber associated 
with that scale. 

In order to solve equations (||) it is standard to de- 
compose physical (perturbed) fields into a spatial and 
temporal part using eigenfunctions which are solutions 
of the Helmholtz equation ||l6| . In the case of the shear 
tensor we write 
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where is a tensor harmonic satisfying 



f( k ) — 3, . n bcd 4- L e bcd r>( k ) a. 



and 
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It is straight forward to verify that the spatial functions 
£? fc ) and also satisfy the Helmholtz equation. 

In order to estimate the dynamical importance of our 
fields we introduce expansion normalised variables 
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giving us a set of scale invariant functions (see, e.g. ||l7| ). 
Introducing the conformal time parameter rj (whose 
defining equation is r] — a -1 ), the scale factor and Hub- 
ble parameter for a dust FRW background are given by 
a{rf) = r/ 2 , H = 2ry~ 3 (see, e.g. ||lq|). 



Substituting these into Eqs. (U2J) we obtain 
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We can also define higher order harmonics by taking co- 
moving spatial derivatives of the lower order harmonics, 
for example Q^ h \ = aV a Q^ can easily be shown to sat- 
isfy 



a 2 ^abc 
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Using the above solution (0) (dropping the index n 
which indicates the scale length of Bft), the decomposi- 
tion (||) and writing 



the wave equations (0) become 
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Equations ( |15[ ) can be solved exactly in the long wave 
length gravitational wave limit (i.e. the term k 2 /a 2 is 
small compared to the other terms in the wave equa- 
tions) , but a numerical investigation gives more transpar- 
ent results. It turns out to be convenient to introduce the 
variables A^ n) = £( fc) rj§/£(?7o), 5(fc,n) = 7Y( fc )?7 2 /E(?7o), 
where rjo and £(r/o) are respectively the initial values of 
the conformal time and the normalised shear, and we 
have reinstated the index n indicating the scale length of 
B§ . These variables are invariant with respect to changes 
in rjo and £(?7o), thus giving us a scale invariant measure 
of the generated electromagnetic field. Performing the 
integration for various values of the wave number n for 
the magnetic test field, we find that the normalised elec- 
tric field A(k tV \ tends to a constant value which depends 
linearly on the initial value of the shear perturbation, 
while the normalised magnetic field B^, n ) tends asymp- 
totically to zero (see Figs. 1 & 2 below). 

In the long wave length gravitational wave case, the ex- 
pansion normalised shear is given by £(fc 1<mg ) = \C\rj 2 — 
3C 2 ?/~ 3 , where C\ and C 2 are integration constants. For 
late times, during the matter dominated era (when the 
equation of state p = applies), the second term in 
£(fc lons ) can be neglected. In this way we can easily ob- 
tain the late time behaviour of the expansion normalised 
gravitationally induced electric field: 
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where A# = is the Hubble radius during that epoch. 
It follows the generated electric field is proportional to 
the expansion normalised shear. Since the normalised 
magnetic field asymptotically tends to zero, the above 
results demonstrate that electric fields produced by this 
effect could play an important dynamical role in the early 
universe, possibly causing charge separation. Further- 
more, because the asymptotic value of the electric field is 
proportional to the magnitude of the shear, we can use 
the CMB anisotropy limits on S to give an upper bound 
on the size of this effect [Eal: 
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VI. DISCUSSION 

In this paper we derived a set of non - linear wave equa- 
tions which demonstrate how electromagnetic fields can 
be generated when gravitational waves couple to an near - 
homogeneous magnetic test field propagating on a FRW 
background. In particular we found that for long wave- 
length gravitational waves, the gravitationally induced 
fields are proportional to the magnitude of the expan- 
sion normalised shear which characterise tensor pertur- 
bations. This allows a simple determination of an up- 
per bound on the magnitude of these fields based on the 
quadrapole anisotropy of the Cosmic Microwave Back- 
ground. 

We note that this paper has not considered the back - 
reaction of this effect on the gravitational dynamics, 
which although small may also give rise to interesting 
results. This issue will be explored in a forthcoming pa- 
per @. 
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FIG. 1. Expansion normalised gravitationally induced elec- 
tric field for different values of the magnetic wavenumber n. 
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FIG. 2. Expansion normalised gravitationally induced 
magnetic field for different values of the magnetic wavenum- 
ber n. 
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